Interplay of spin-orbit coupling and Zeeman splitting in the absorption lineshape of 
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We suggest that electron spin resonance (ESR) experiment can be used as a probe of spinon exci- 
tations of hypothetical spin-liquid state of frustrated antiferromagnet in the presence of asymmetric 
Dzyaloshinskii-Moriya (DM) interaction. We describe assumptions under which the ESR response 
is reduced to the response of 2D electron gas with Rashba spin-orbit coupling. Unlike previous 
treatments, the spin-orbit coupling, Aso, is not assumed small compared to the Zeeman splitting, 
Az- We demonstrate that ESR response diverges at the edges of the absorption spectrum for ac 
magnetic field perpendicular to the static field. At the compensation point, Aso ~ Az, the broad 
absorption spectrum exhibits features that evolve with temperature, T, even when T is comparable 
to the Fermi energy. 
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I. INTRODUCTION 

In a spin system with isotropic exchange the absorp- 
tion spectrum of ac magnetic field is a S-peak at w = Az, 
where Az is the Zeeman splitting, independently of the 
exchange interaction strength^. Similarly to the Kohn 
theorem for cyclotron resonance^, this fact is the con- 
sequence of coupling of spatially homogeneous exciting 
field to the center of mass of the system. Therefore, any 
deviation of the absorption spectrum from the ^-shape 
implies a violation of the spin-rotation symmetry which 
is either due to anisotropic terms in the Hamiltonian or 
due to development of spontaneous ordering below criti- 
cal temperature. 

The subject of the present paper is the shape of the ab- 
sorption in the spin-liquid ground state'^ of the frustrated 
antiferromagnet in the presence of Dzyaloshinskii-Moriya 
(DM) interaction which originates from the spin-orbit in- 
teraction in the underlying system of electrons^i^. We 
find that electron spin resonance (ESR), which measures 
the ac absorption, can be used to probe the fractional- 
ized spinon excitations. Specifically, we assume that the 
ground state of the frustrated spin system is described by 
the Fermi sea of the neutral spin excitations (spinous). 
Emergence of the Fermi sea of neutral excitations is at 
the core of the RVB state proposed by Anderson 25 years 
agc^'-. Recently the spin liquid state with spinon Fermi 
surface was found for half-filled Hubbard model on tri- 
angular lattice^ii^ for intermediate U/t values. 

Dzyaloshinskii-Moriya interaction that gives a nontriv- 
ial shape to the spin resonance is most conveniently intro- 
duced into the Hubbard model via spin-dependent hop- 
ping term 

H = ^{ct^ai^t'^aP + iK] ■ Safs)CjJ3 + h.C.} 

+U^ni^ni^. (1) 

i 

Here 2s is the vector of Pauli matrices, a and /3 are spin 
indices, and \ij is the DM vector on the link [ij). The 
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FIG. 1: Energy dispersion of two subbands in the presence 
of Zeeman splitting and spin-orbit coupling are plotted from 
Eq. ^. The curves cross at k^ = Az/2a. Yellow regions 
designate the momenta of the states taking part in absorption. 
For Fermi-level position Ep = A?z/^i^\, where u>i = 2mo? , 
there is no threshold for absorption. 



first (t) and the last [U) terms are responsible for the 
Fermi surface formation. Here we argue that, in course of 
reduction of these terms to the Hamiltonian of spinous, 
as outlined below, the second term in Eq. ([1]) gener- 
ates a term in the spinon Hamiltonian which has a form 
of spin-orbit coupling familiar from 2D electron gas in 
semiconductor structuresA^.. 

In order to illuminate the central idea of our work 
and to keep technical complications to a minimum, we 
will focus on the case of a uniform spin-orbit interac- 

,). Here vectors r^j denote 



tion Aij- — Xfi X {fi 

nearest-neighbor lattice sites and n = z is a normal to 

the plane. Note in passing that, for a square lattice. 



precisely this arrangement is known to realize in, e.g., 
YBa2Cu306+a~^"— . For the chosen form of spin-orbit 
interaction the DM term in Eq. ([T]) assumes the form 



HsoOi) = -2A ^ cl^^isx sm[ky] - Sy sm[kx]}ck,i3. (2) 

k 

One can recognize in Eq. ([2]) the lattice version of the cel- 
ebrated Rashba spin-orbit term^°'^*' in the Hamiltonian 
of a 2D electron gas. 

Technically, the reduction of the full electronic Hamil- 
tonian, Eqs. (HI) and ©, to that describing spinon exci- 
tations of the uniform spin-liquid stated is achieved with 
the help of the slave-rotor techniqueii"— . Within this 
technique, a spinon fi^a with spin projection a at i-th 
site inherits the spin of the electron, c^.q, while the rotor 
variable, e'^*, inherits its charge, Ci^a — fi,aS~'^^'. Fol- 
lowing the approach initially outlined in Ref. U^ one 
then substitutes this parametrization into the original 
Hamiltonian Eq. ([ij and decouples the resulting cou- 
pled dynamics of spinous and rotors into a sum of two 
separate Hamiltonians describing spinous and rotors in- 
dependently. The parameters of these Hamiltonians are 
determined self-consistently and are renormalized with 
respect to their bare values. In particular, one finds (see 
Sect. HI B 1 of Ref. (fill )) that spinous are described by 
the free fermion Hamiltonian 



Hf = 



iJ 






(3) 



the parameters of which - t°^ , X'^^ - depend on the rotor 
Hamiltonian, Hg, the precise form of which is not needed 
here (it is given by Eq. ([M| of Ref. (fill )). For example. 



tij (cos[i 



„j,j\„„„^„i ffj])^, where the average is taken with 

respect to Hg. Since A^j represents the spin-dependent 
part of the hopping integral, ti_j, we expect it to renor- 
malize similarly to the latter quantity. Note in passing 
that in many cases spin-orbit interaction can be gauged 
away altogether by an appropriate unitary rotation, see 
Ref. [l^ In these cases the above statement becomes 
exact. Note also that the symmetry alone, namely the 
oddness of the spin-orbit interaction in Eq. ([1]) under the 
spatial inversion, dictates that the spin-orbit interaction 
of spinous has the same functional form. 

The last term in Eq. ([3]) represents Zeeman coupling 
of electron's spin to the external magnetic field, h, which 
does not involve charge degrees of freedom at all, and 
thus preserves its form for spinous as well. In gen- 
eral, magnetic field also couples to the orbital motion 
of spinous via the higher order closed loop processes^ 
(ex t^/U'^ for triangular lattice, for example). For the 
sake of argument we do not consider these processes here: 
for example "orbital" coupling is absent in geometry 
where magnetic field is parallel to the two-dimensional 
layer. 

We would like to emphasize that the Hubbard interac- 
tion, U, does not enter Eq. (|3]) at all. In the slave-rotor 



representation it is present only in the rotor Hamiltonian 
Hg and affects spinous indirectly, via renormalized pa- 
rameters such as f^^ in Eq. ^, see Refs. [ll| - [l3l This 
feature of the slave-rotor theory makes it clear that, as 
long as the spin-orbit and Zccman interactions are small 
perturbations to electron's Hamiltonian (which, practi- 
cally speaking, is always the case), the spinon sector of 
the Hubbard model must be described by Eq. ([3]), so 
that the spatial structure of the spinon spin-orbit cou- 
pling, A^*^, inherits that of the electron spin-orbit term 
in Eq. (|1|). Here it means that spin-orbit part of Eq. 
(|31) is given in the momentum space by Eq. ^ when 
electron operators Ck,a are replaced by spinon fk^a ones. 

Finally, the spin-liquid phase we are interested in is 
just a disordered (Mott) phase of Hg, in which the ex- 
pectation value of the rotor field is zero, (e*^) = 0, and 
the charge degrees of freedom are gapped. (The phase 
with a condensate of charge rotor field (e*^) 7^ describes 
the usual metallic phase of Eq. ([T]) and is not of inter- 
est here.) Importantly, spin excitations of this phase are 
neutral spinons, fi^a, forming the Fermi-liquid state and 
described by the Hamiltonian Eq. ([3]) with parameters 
which are renormalized by the gapped charge fiuctua- 
tions (the notion that spinons can be treated as neutral 
fermions^iiiii^ has an important caveat that they inter- 
act strongly with emerging gauge fiel d^'^'^^ ). 

We thus see that, under the assumptions described 
above, the Hamiltonian of strongly interacting Fermi sys- 
tem Eq. (|T]) in the spin-liquid phase maps onto that of 
non-interacting spinon Fermi gas with spin-orbit interac- 
tion of Rashba type ([3]). Hence, the basic features of the 
ESR response in an exotic spin liquid state can be cap- 
tured within a much simpler model of two-dimensional 
electron gas subject to spin-orbit interaction. Surpris- 
ingly, despite several decades of intensive research in the 
latter area, we are not aware of the calculation of ESR 
response in two-dimensional electron gas subject to both 
spin-orbit and Zeeman interactions. 

In earlier studies of spin-orbit coupled two-dimensional 
electron gas the combined effect of Zeeman field and spin- 
orbit interaction was considered in relation to electric 
dipole spin resonance^S, (EDSR), which is the absorption 
of the ac electric field between the Zeeman-split levels. 
Obviously, it is spin-orbit interaction which allows this 
absorption to occur, since it couples the electron spin 
to the electric field. In theoretical papers^i^ on EDSR 
in two-dimensional electron gas a spin-orbit term of the 
form Eq. 1^ in the small-k limit was added to the Hamil- 
tonian of free electrons. Then the matrix element of tran- 
sition between two spin levels was calculated for differ- 
ent orientations of magnetic field. However, we cannot 
use the results of Refs. [23ll24l for finding the lineshape 
of ESR. This is because in these papers it was assumed 
that the magnitude, A50, of the spin-orbit term is much 
smaller than the Zeeman splitting, A^, so that Aso was 
treated perturbatively. As a result, in the absence of 
disorder, the lineshape of the resonance was simply a S- 
peak. The natural basis in which the magnitude of this 
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FIG. 2: (a) Two Fermi surfaces corresponding to Ep — 
/S.z/4ijj\ and oji/Az = 0.1 are plotted from Eq. (|30p . where 
the projections of momenta kx and ky are in the units l^z/'2ct. 
(b) Evolution of the Fermi surfaces near the point ky = Q 
for different Ep (enclosed by the gray box in (a)): (thin 
line) Ef = 0.95S^°', (thick hue) Ef = E^ , (dashed hue) 
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peak is calculated2^>2i is the basis of spin projections on 
the applied magnetic field. 

Another group of relevant earlier paper a^^'^^ dealt with 
the phenomenon of chiral resonance, which takes place 
when A50 is finite. In this case spin-orbit term splits 
the free-electron spectrum into two branches, separated 
by Ago J which correspond to different chiralities. Then 
the dipole absorption of the ac electric field between the 
branches is allowed. The absorption spectrum of the chi- 
ral resonance has a box-like shape^ - centered at w = A50 
with a width of the order of A'^q/Ef, where Ep is the 
Fermi energy. The absorption is calculated in a natural 
basis of spinors, describing opposite chiralities. Still, we 
cannot use the results of Refs. .25„2& for calculating the 
ESR lineshape, since magnetic field was assumed zero 
in these papers. Instead, we will adopt the general ap- 
proach of Ref. i26| to the calculation of the ac response. 
This approach is based on the calculation of the opti- 
cal conductivity. We will subsequently demonstrate that 
the ESR lineshape can be expressed through the optical 
conductivity in a simple way. 

The outline of the paper is as follows. Calculation of 
the optical conductivity for arbitrary relation between 
A50 and Az, is presented in Sect. II of the present 
paper. In Sect. III. the results of this calculation are 



used to find the ESR lineshape for different orientations 
of the ac magnetic field. 

Normally, (t^lj) is a peak with two sharp edges, which 
are imposed by the energy conservation, as illustrated 
in Fig. [TJ Notably, we demonstrate that (7{uj) exhibits 
a singular behavior near these edges. The character of 
singularities is different for ac field parallel and perpen- 
dicular to the external in-plane magnetic field. We trace 
how these singularities are smeared out upon increasing 
the temperature, T. 

The most nontrivial behavior of cr(a;) corresponds to 
the situation Az = Ago, when Zeeman splitting is "com- 
pensated" by the spin-orbit splitting along the direc- 
tion of the magnetic field. As a result of this compen- 
sation, the two branches of electron spectrum, £(^'(k) 
and £*^^^(k), cross each other at certain momentum k — 
{ko,0) while the Fermi level Ep = S^^^k) = £:(2)(k) is 
located at the point of crossing, as illustrated in Fig. [T] 
The underlying reason why the shape of o'(uj) is peculiar 
near the condition Az = Aso is the following. 

As the difference Az — Aso changes sign, the Fermi 
surfaces ^^^•'(k) — Ep and £'^^)(k) — Ep undergo rapid 
restructuring. This restructuring is illustrated in Fig. [21 
Since the states responsible for absorption lie between 
the two Fermi points, see Fig. [U in the vicinity of Az = 
A50 the absorption remains finite at low frequencies to <C 
Az- We will show that in this regime the absorption 
spectrum exhibits features that evolve with temperature, 
even when T is comparable to the Fermi energy. 



II. OPTICAL CONDUCTIVITY WITH FINITE 

Az AND Aso 

A. Hamiltonian and eigenfunctions 

Assume that electrons are confined to the {x, y) plane, 
while magnetic field, creating the splitting, Az, is di- 
rected along the y-axis. For concreteness we choose the 
Rashba-typei2i2£ form of spin-orbit coupling 



i/so(k) = 2a{k, 



XJy 



i^y^x )■> 



(4) 



where a is the spin-orbit constant. In the matrix form, 
the Hamiltonian of the system reads 



2m (^x 



i-s2- — a{ikx + ky) 



H = 



V' 
"*~2 Ol.\~1'kx + ky) 2m^ -c ~^ y' 



,(5) 



where m is the effective mass (we set fi = 1 through- 
out the paper). The two branches of the spectrum and 
corresponding spinors are given by 






^1(2) 



afc^ - ^ ) : (6) 



u 



(1),(2) 



1 I fc^ c-i(2 



V2 



(7) 



V 



1 



To express the optical conductivity in terms of the spec- 
trum Eq. (|6|) and eigenfunctions Eq. Q one has to add 
to the Hamihonian Eq. ^ the term 



Hi, 



eE • p cos Lot, 



(8) 



where Ecoswi is the external ac electric field, and cal- 
culate the ac current. This calculation is not straightfor- 
ward because the eigenfunctions from different branches 
corresponding to the same momentum are orthogonal to 
each other. For this reason, it is much more convenient 
to adopt the approach of Ref. [23. Within this approach 
the system response to the scalar perturbation 



Vip, t) = Vo cos{iot - q • p). 



(9) 



with finite wave vector, q, is studied. The energy absorp- 
tion rate, /(q), caused by this perturbation, is calculated 
using the Golden rule. The optical conductivity at finite 
q is related to /(q) as follows 



cr(q,a;) 



1 



^/^)/(q) 



2Vo'q 



2„2 



(10) 



where exp(— w/T) accounts for the emission processes. 
The static limit q ^ is taken only at the final step 
yielding 

.H..e^(l-e-W-)lim(^)^|(.W|4^)^)|^ 



'.S{£, 



(2) 
k+q 



4^^-: 



i-/(4+J/(^i 



^l^') 



(11) 



where f{£) = (l + exp[(f - Ef)/T]\ is the Fermi 
distribution. 

The formal reason why absorption of the ac field is 
permitted in the presence of spin-orbit coupling is that, 
at small q, the product (u^ l'"k+ 
of q, as readily follows from Eq. (O 



u^ ' ) is a linear function 



i'VX) 



./a\ {ak^ - ^Az)qy - akyq^ 



\2) {ak, - \AzY + a-'kl ' 



so that the combination, | (u^. 



^^^^^ 



^/(?^, remains finite 



in the limit of q — > 0. In subsequent sections we analyze 
Eq. ([TT|) first at T = 0, and then at finite temperatures. 



B. Optical conductivity at T = 

It follows from Eq. p^ that at finite Zeeman splitting 
the optical conductivity is highly anisotropic. For polar- 
ization perpendicular to the magnetic field, vector q is 



directed along the x-axis. Substituting Eq. p2|) into Eq. 
(|lip. and converting the sum into an integral in polar 
coordinates, we get 



crx(w): 



xe 



2 4 

e a w 



167r 
fc2 



n{4>Y 



2m 



dkt 



n{(t)) -Ef]Q\Ef 



sin^ (j) 5{2n{(f)) - Lo) 



2m 



where we have introduced the notation 
17(0) — \ a^k"^ — a/S.zkcos(j) - 






(13) 
(14) 



The factor £7(0) ** in the integrand emerges from the de- 
nominator in the product {u 



Wi 



k I "k+q. 



, seeEq. (HH). The 
numerator in this product is proportional to fcj,, giving 
rise to sin in the integrand. It is convenient to first 
perform the angular integration in Eq. (J13p with the 
help of the (5-function. Then the remaining integral over 
k can be cast in the form 



a^{uj) = 



/tFo; A^ 
X I dkk 



y/m{2Ef 



a^k^Al - a^k^ 



A^ 






1/2 



(15) 



Integration in Eq. (jisp can be performed explicitly, yield- 



ing 



(7^{uj)-- 



Svr 



dry/l 



—— I ( arcsinr -I- r\J\ 
Ibn , 



(16) 



The limits of integration in Eq. (J16l) are given by 



rfA^) 



1 



:(Al 



2Azc.^^^o 
where the frequency cji is defined as 

A^ 
AEf 



a;2±2wiw-A|), (17) 



(18) 



Expression Eq. (1161) is valid when both r^, rf lie within 
interval —1 < ri,rf < 1. For frequencies at which both 
ri and rf are outside this interval, we have a±{uj) = 0. 
When only one of the limits, say rf, is outside the inter- 
val, one should set rf = 1 in Eq. p6|) . 

For polarization parallel to the magnetic field, the cal- 
culation of conductivity is completely similar to the cal- 
culation of (T+(a;). The only difference is in the product. 



(■"k l^k+q)' which, for q directed along y, is proportional 
to {akx — Az/2). In polar coordinates this factor leads 
to replacement of sin^ in the integrand of Eq. (|13p by 
(cos0 — Az/ak)^. As a result, instead of Eq. (fT6)) . the 
momentum integration reduces to 



aJuj)-- 



\ frf 



dr- 



2 

' arcsm r — r 



Wit I V 



xA 



(19) 
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Equations p^ and p^ constitute our main results. We 
start their analysis from the limiting cases Aso ^ Az 
and Aso < Az- 



1. Aso:^Az 

It is instructive to trace how the shape of the chiral 
resonance emerges from Eqs. ([T6|) . ([T9l) . Taking the limit 
A^ -^ 0, we find that r^ and r/ should be set to —1 and 
1, respectively, for w within the interval |a; — Aso\ < wi. 
Outside of this interval we have ri = rf — ±1 and the 
absorption is zero. This yields the lineshape 



cr(a;): 



16/ 



8 



A2 



2{jjuj\ 



e 



A2 



1UJUJ\ 



(20) 



with sharp ends, which can be further simplified by in- 
troducing bix) = bj — Aso and taking into account that 
wi <C Aso- Then we reproduce the result of Refs. l26l)28l 



16 \ wi 



(21) 



where Q(x) is the step-function. 

Next we study how a small Zeeman splitting, Az <C 
wi, smears the steps at w — Aso = ±^i- Introducing 
dimensionless deviation from the right boundary 



(A 



so 



uj -t- Wi) 



?7(w) 

we get the following behavior of uj 
boundary 



(22) 



(Til near the right 



167r 



- - arcsin77(i:j) =F v{^) V^~rf{uj) 



(23) 
This behavior is illustrated in Fig. [3l We see that a step 
at 5lj — wi extends, at finite Az, over an interval 6ui — 
wi = ±Az. Optical conductivity is strongly anisotropic 
within this interval. In particular, an has a zero slope at 
the center of the interval. As Az exceeds oji, the flat top 
in optical conductivity disappears completely. 



FIG. 3: Smearing of the right edge of the chiral resonance 
at small Az <C t^i is plotted from Eq. (|23|) as a function of 
dimensionless deviation, r]{uj), Eq. (|22ll . from the boundary. 
(j|| and crj_ are plotted in the units e^/16. 



2. Az > Aso 

Consider now the EDSR limit of strong Zeeman split- 
ting, Az > Aso- Then Eqs. ^, ^ reduce to 



o"±(w) = 



iM 



/V2e2\ 




\Alo~{u:-Azf 


1 4-; 


2Azuj 


^' ^ 




A|o-(a.-Az)21 


4V^7ry 


2A2W 



-,1/2 



(24) 



(25) 

The fact that the EDSR spectrum is confined to the in- 
terval |w — Az\ < Aso could be expected on general 
grounds. It is also natural that optical conductivity is 
proportional, via wi, to the square of the SO coupling 
strength. Less trivial is the fact that EDSR absorption 
is highly anisotropic, so that a± turns to zero at the 
edges, while cr|| (uj) exhibits a divergence near the edges. 
This can be seen in Fig. |3^. 



A. 



A 



so 



Consider now the intermediate regime when Az and 
Aso are of the same order. Since we assumed that Fermi 
energy is bigger than both Az and Asoj the condition 
Az '^ Aso automatically implies that both Az and Aso 
are much bigger than uJi. Then the integrands in Eqs. 
(fTHll . (fT9|) can simply be multiplied by rf — ri = 2uJi/Az, 
yielding 



aj_{uj) 



2uji 



^JA^^ 



r'^iu), 



-iiH=(^ 



2wi r2(cj) 
A^^l-r2(^' 



(26) 



(27) 
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by a maximum in a±. For opposite relation, A^ > A50, 
there is a minimum in r^(a;) at w = cJq = (A^ — A|q)^/^. 
This minimum translates into a minimum in a\\ {ui) and a 
maximum in cr^(a;). The behaviors of cr\\{ijj) and cr_L(w) 
for different relations between A^ and A50 are illus- 
trated in Fig. 131 

The most interesting behavior of optical conductivity 
takes place near the point of "compensation" , where A^ 
and A50 are close to each other. Then the minimum 
in an and the maximum in a± become very sharp. This 
implies that absorption becomes possible at very low fre- 
quencies. On the physical level, the low-frequency ab- 
sorption is the result of crossing of the Fermi surfaces 
£:(i)(k) = Ef and £(2)(k) = Ep, see Fig. [H which takes 
place at A^ = Ago- The Fermi energy is equal to 



E 



(0) 



A2 



(29) 



at the point of crossing. As the difference Az — Ago 
increases the Fermi surfaces undergo restructuring. In 
Fig. [5] we show the evolution of the Fermi surfaces 



^- - ^ (^^ 



fc,1 ± 



I oi^k'^ 



akx 



A^ 
2 



(30) 
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as Ep departs from Ep ' . Since the states responsible for 
absorption lie between the two Fermi points. Fig. [U the 
low-frequency part of the absorption spectrum follows 
this departure. 



C. Optical conductivity at finite temperature 

We now turn to the general expression Eq. pT|) for 
optical conductivity and perform integration over the di- 
rections of momentum, k. This yields the following gen- 
eralizations of Eqs. ([TB)) . p^ to finite T 



FIG. 4: Optical conductivity at zero temperature, Eqs. (|26|) . 
(|27|l vs. dimensionless frequency, uj/Az, are plotted in the 
units (e^ /it){uji/Az), for three values of the ratio Aso/Az- 
(a) Aso/Az = 0.5, Zeeman splitting dominates, (b) Zeeman 
and SO splitting nearly "compensate" each other, Aso/Az ~ 
0.99. (c) Aso/Az = 10, SO splitting dominates. 



a I (oj) = = — : — 1 sinh , 



(^\ 



dk k 



Aa'^k^ ^ [lo + AzY {uj - AzY - Aa'^k 



cosh I 



2mT 



:)+cosh(#) 



(31) 



where 



r(w) = 



2Azuj 



(A|o - -^ - Al 



(28) 



is given by Eq. ((T7| with uji — Q. It is seen that 
CTx monotonically grows with r^, while a\\ monotoni- 
cally decreases with r^. At the same time, the fre- 
quency dependencies of ctx a-nd crii are strongly non- 
monotonic. Firstly, for A50 > A^, r{uj) turns to zero at 



UJ 



Wo 



(A|o 



A^)-'^/^. A zero in a\\ is accompanied 



cr||(w) 



— T= — - — i sinh ( — - ) / dk k 
2^ujAz) \2T) Jo 

(4a2fc2 - A| 



4a2fc2 - (tj + Az)2 (a; - Azf - ia^k^ 
1 



cosh (513.-^)+ cosh (i^) 



(32) 
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we find that, outside the plateau, (j{(jj) falls off as 
exp(— ^'^ \5uj ± Wil). Thus, the plateau is only slightly 
rounded near the edges, namely, within the interval 
\6uj ± wi| ~ Tuji/Aso ^ ^1- Upon increasing T, at 
T ^ Ago, we can neglect 1 in both cosh-terms. Then 
the denominator will become cosh {Aso<>i-^/4:Tuji). We 
see that characteristic width of the absorption spectrum 
is ~ TuJi/Asoj which is much bigger than wi. Note 
however, that, since wi = A^q/AEp, this width remains 
smaller than Ago as long as T remains smaller than Ep. 
Evolution of the shape of the chiral resonance with T is 
illustrated in Fig. [S] 



2. Temperature smearing of EDSR 



FIG. 5: The lineshape of the chiral resonance is plotted from 
Eq. (|33p in the units of e^/16 for three dimensionless temper- 
atures, Aso/T, and spin-orbit coupling strength Aso/Ef — 
0.01. In the low-temperature limit, Aso/T = 50 (thick line), 
the lineshape is box- like. At temperature Aso/T = 15 (thin 
line) the edges of the box get smeared. Complete smearing 
corresponding to Aso/T = 2 is illustrated with dashed line. 



Expressions, Eqs. (|31l) and ([5^ . describe the broadening 
of the optical conductivity peak at finite temperature. 
We can expect on general grounds that the chiral reso- 
nance broadens with T much faster than the EDSR peak. 
This is because the distance between the two Zeeman- 
split branches of the spectrum is Az at any momentum, 
while the SO-splitting is proportional to the momentum. 
In the latter case, the transitions between the states away 
from the Fermi surface, which become possible with in- 
creasing T, have different frequencies. For the same rea- 
son, the low-frequency absorption that becomes possible 
near the condition Aso = ^z is most sensitive to in- 
creasing T. 



1. Temperature smearing of the chiral resonance 

In particular case of the chiral resonance, Az -^ 0, the 
integration over k can be performed explicitly, since the 
integrand is non-zero only in the vicinity of A: = uj /la. 
Then it is sufficient to substitute fc = uj /la into the de- 
nominator. This leads^^ to the following finite-T gener- 
alization of Eq. (|2T]) 



-(MH^^ 




At low T <ti Aso the hyperbolic sine and cosines can 
be replaced by the corresponding exponents. Then 



Consider now the opposite limit of EDSR. The prime 
modification which finite but low temperature imposes 
on T = expression Eq. p4)) is a smearing of the 
square-root singularities at the edges. Due to temper- 
ature smearing the factor (|w — Az\ ± Aso)^^'^ gets re- 

1 /2 

placed by {uiiT/Aso) ■ To see this, it is convenient to 
transform from integration over k to integration over r, 
as it was done at T = 0. Then from Eq. pTjl we get 



CT± 



dry 1 — r'^ 
sinh(a;/2T) 






j + cosh(a;/2T) 



.(34) 



It is seen that for T <C Az we reproduce the zero- 
temperature result. In the opposite limit, T ^ Az, the 
major contribution to the integral still comes from r close 
to —1. In this limit, the hyperbolic sine and cosine can be 
replaced by the argument and 1, respectively. Introduc- 
ing the dimensionless variables, v = {Azuj/2ujiT){r -I- 1) 
and 



a{ijj) = 



[Az^ujf 



A2 



2Azuj 



(35) 



we can rewrite tr i as 



-M^y £ 



3/2 



T 
A^ 



1/2 



5J 



|a(u;)|A^u; \ 
2uJiT ) ' 
(36) 
where the dimensionless function g± (x) is defined as 



9i^{x) 



dv^/v 



cosh(a; ±v) + \ 



(37) 



In the above definition the sign is minus for a{uj) > 
and plus for a(w) < 0. Consider first positive a(w). For 
X ^ 1, the function g_L{x) behaves as 2x^/^, and we 
reproduce the zero-temperature result. For x ^ 1, we 
can replace g^{x) by ^(0) w 1.07, so that 



"-^='4^ 



Az 



3/2 



T 

A^ 



1/2 



(38) 



Comparing Eqs. (1371) and (p4|) . we conclude that the 
smearing of the square-root singularity takes place within 
the interval, (A^ — a; — A50) ^ ^iT/Aso- 

The large-x asymptote, for negative a, of g± is gj, = 
•\/7re~^. The exponential describes the decay of a± in the 
region of lo where it was zero at T = 0. 

In a similar way, the smearing of inverse square-root 
divergence in Eq. (P5|) is described by the function ^n (x) 
defined as 



9\\{x) 



Then Eq. ([32]) reads 



dv 



a«{uj) = 



y/v(^ cosh(x ± u) -I- 1) 



= f^\ /cJi\i/2 f \a{Lj)\Azuj 



\16tt 



T 



2ujiT 



(39) 



(40) 



For positive a(w) in the limit a; S> 1, the function gn 
tends to 2/^ and Eq. (gO]) reduces to Eq. ([25]). The 
inverse square-root singularity is smeared in the same 
domain, (A^ — w — A50) ~ ujiT/Agoi as for cr^, where 
the argument of g\\ is '^ 1. In this domain the ratio a\\ /a± 

remains big but finite, aii/a± ^ — S^. For a{uj) negative, 
in the large- a: limit, we have gi\ = 2ypKe~''^ . Substituting 
this ^11 into Eq. (HUl) . we find that to the left of the 
smeared anomaly an decays as 



Uv^y ^T 



1/2 



exp 



^A 



so 



2uJiT 



Aso - (w - A2 



(41) 



3. Az ~ Aso cit finite T 



As stated in the Introduction, the most nontrivaial 
situation corresponds to Az = Aso, where the two 
Fermi surfaces undergo restructuring. The expressions 
for T — optical conductivity under this condition read 
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(42) 



(43) 



We will now analyze how temperature affects the low-w 
part of absorption and the sharp edge at a; = 2A^. 

It is convenient to make the change of variable A:^ = 
2m{EF + QT) in Eq. dST]), after which it takes the form 



cr±(cj) = 



327ra;2 
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so 
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FIG. 6: Evolution of optical conductivity with T in the 
regime Az = Aso- Eqs. 1)45^ and H46p . are plotted in the 
units {e^ /1Q-k)(Az / Ef) for three values of T/Ef\ (a) low 
temperature limit T/Ef = 0.05. (b) intermediate tempera- 
ture T/Ef = 0.2. (c) "high" temperature T/Ef = 0.5. Inset 
in (a) shows the dimensionless functions G^{a) and G\\(s) 
calculated from Eqs. (|49ll and (|50[) . respectively. 



Now look at the point where Aso ~ Az, in the limit 
T > w. Then Eq. gH) reduces to 



o-_l(w) 



(44) 



e 
64^ 



luEf 



dQ 
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(QT 
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cosh Q + 1 



(45) 



Similarly, the T ^ oj limit of Eq. (|32|) assumes the form 



crii (w) 



f. 



\64:TT J \ijEp 
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Consider first the domain uj <C A^. Then in the com- 
bination (^ ^) in the integrands of Eqs. (HSl) and 

(|^S)) . the term (w/A^)^ can be dropped. This allows us 
to simplify Eqs. (gS]) and (05]) to 



a±{uj)-- 



e 
32^ 



Az 
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32'7r 



^l|H=(:^)(;^ 



/2w£;f\ 



AzTV ' 

where the functions G± and Gh are defined as 
1 f" dQ 



G^is) 



G||(s) = 



1 



_g coshQ + 1 



dQ 



Vs 



Q' 



(47) 



(48) 



(49) 



(50) 



L, coshQ + 1 ^s2 _ g2 

These functions are plotted in Fig. [S^ inset. At s ^ 1 
they both behave as 7rs/4. For s :^ 1 the function G± 
saturates at G_l(oo) = 2, while the function G\\ falls off 
asG||(s)«27rV3s2. 

It can be seen from Fig. [5^ inset that G||(s) has a 
maximum at s = s = 1.83. This translates into the 
maximum in absorption at frequency 

. O.QlAzT 

This sharp maximum clearly shows in the absorption 
lines plotted in Fig. [S^,b which were calculated numer- 
ically from the full expression Eq. (I46p for particular 
temperatures, T/Ep = 0.05 and T/Ep = 0.2. Remark- 
ably, the position of maximum moves linearly with tem- 
perature for low temperatures. The maximum starts 
to disappear when temperature becomes high enough, 
T/Ef = 0.5. 

It is also seen from Figs. [S]d,c that the right boundary 
oj sa 2Az of the absorption spectrum is smeared with 
temperature. This smearing can be described analyti- 
cally in terms of the functions g± and 511 defined by Eqs. 
([57)1 and ([5^ . Indeed, for uj close to 2Az, the square 
root in the numerator of Eq. (P5|) can be simplified to 



2 

Az 



[UJ 



V2iu;/Az)[^ 

aj_{uj) acquire the form 



cr±{uj) 



2A; 



,1/2 



Then crii (a;) and 




cr||(w) 



A; 



\64:V2t:J V Ef 



El 

T 



1/2 



2Ef {lo ~ 2Az] 



T 



A, 



(53) 
It is worth noting that even at high temperature T — 
0.5Ef the right boundary uj = 2Az is very well pro- 
nounced. 



III. ESR LINESHAPE 

In this section we use the above results for the absorp- 
tion of the electric field to calculate the ESR lineshape. 

For ESR, the corresponding interaction Hamiltonian 
of spins with the ac magnetic field reads 



H^, = 7(2s • B) cos ujt, 



(54) 



where 7 is the gyromagnetic ratio. The quantity describ- 
ing the ESR strength, imaginary part of susceptibility, 
can be expressed via the eigenfunctions of the free Hamil- 
tonian Eq. ([S|) as follows 



Xi(w) = 87r7^(l -e" 



k 



,(1) 



k |*j|'"k 



(2) 



x6 £, 



(2) 



7(1) 



^)[l-/(^l 



(2)^ 
k , 



/(^k 



(1)^ 



(55) 



where i is the polarization axis of the ac magnetic field. 
Our prime observation that the values Xi can be ex- 
pressed in a simple way though the components of op- 
tical conductivity, a±{uj) and cr||(w), calculated above. 
Consider, for example, ac magnetic field polarized along 
X. Then, using Eqs. (jS]), ©, the matrix element 



,(2)|s L,(i)\ 



can be expressed as 



/ (1)1 , 



(2)\ _ . gfcj; - Az/2 



£. 



(2) 



£ 



(1) 



(56) 



Upon transformation to polar coordinates, the numerator 
of Eq. ([56| turns into k cos </)— A2/2, while the (5-function 
in Eq. (j55l) insures that the denominator is equal to uj. 
In this way the summation over momenta in Eq. (j55p 
reduces to same integration as in Eq. ()19|) and we arrive 
at the relation 



r ^ 27^ , . 

Xx[^) = -i^-^uja\\{uj). 



(57) 



From comparing Eq. ([55t to Eq. ([TTj). it is obvious that 
this relation applies at all temperatures. 



Similarly, the matrix element 



,(1) 



Sv\U 



(2) 



) in polar 



coordinates assumes the form iak sin (l)/{£\ 
makes Eq. ([55)1 for 
and leads to 



-£. 



(1)^ 



k '-'k ^ 

y proportional to a±_ in Eq. 



This 



Xvi^ 



^2^--^(-)- 



(58) 
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Concerning polarization of magnetic field along 2, we 
have 



(4^^ I 



,(2)\ _ 



-1, 



(59) 



so that summation over momenta does not contain an- 
gular dependence at all. This allows to relate Xz{^) to 
the sum of (T_i_(uj) and cr\\{(^) as follows 



Xz = -^-^^ [c^± {^) + c^ll i^)] ■ 



(60) 



We see that ESR line shape for different polarizations of 
the ac magnetic field can be written in terms of optical 
conductivity a± and crii . Response to the ac field po- 
larized perpendicular to the static field, namely, Xx and 
Xz, is determined mostly by cr||(a;) which diverges in an 
inverse square-root fashion at the boundaries of the ab- 
sorption interval. In weak static field the lineshape is a 
(5-peak at w = Aso^- 

It is natural that the ac field parallel to the static field 
causes less singular response, than the ac field normal 
to the static field. Indeed, Xy is proportional to a±(uj), 
which exhibits square-root suppression at the boundaries 
of the absorption interval. 



IV. DISCUSSION 

Interplay of spin-orbit interaction and in-plane mag- 
netic field results in strongly anisotropic line shape of 
ESR response of the two-dimensional fermion gas. The 
response is most singular at the upper and lower bound- 
aries of fermion continuum when the ac magnetic field 
is perpendicular to the static magnetic field, see Fig. |4] 
and Section IIIIl The behavior of the ESR lineshape is 
most "lively" near the "compensation" point Az — Aso, 
where it exhibits a nontrivial temperature dependence up 
to high temperature T ~ Ep. 

A possible way to reveal these features experimentally 
is to fix the polarization of microwave radiation and to 
rotate the sample with respect to the external static mag- 
netic field. 

Restricting two-dimensional gas to a quantum wire ge- 
ometry when the motion along, say, y-axis is quantized, 
will result in "slicing" the Fermi-surface along several al- 
lowed ky momenta^-. In the extreme case of the single 
channel quantum wire, when only the lowest transverse 
subband is populated, the ESR lineshape will be reduced 
to two delta-function like features at the lower and upper 
boundaries of the ESR continuum. This is exactly the re- 
sponse found previously^ within a different, purely one- 
dimensional treatment of an insulating Heisenberg spin 
chain with uniform DM interaction. Ref. 32 finds a very 
similar response of a conducting quantum wire with spin- 
orbit interaction to ac electric field. The correspondence 



between the two above mentioned one-dimensional stud- 
ies illustrates nicely the central point of our work that 
ESR is an effective probe of the Fermi surfaces irrespec- 
tive of their origin. Both conventional Fermi surface^ 
and spinon Fermi surface^ give rise to the same ESR 
response. 

One-dimensional geometry also allows one to address 
effects of electron-electron interactions omitted in the 
present work. Although relatively marginal even in one- 
dimensional geometry^ij these are difficult to account 
for at finite frequencies as recent detailed investigation 
shows^. Extending this approach to two-dimensional 
systems, along the lines of renormalization group calcu- 
lations in^, represents an interesting outstanding prob- 
lem. 

Throughout the paper we neglected the disorder. In re- 
alistic sample, with finite disorder, the absorption peak 
calculated above will always reside on the background of 
the Drude peak coming from the intraband scattering. 
The condition that disorder does not smear the inter- 
band absorption peak reads: A^t, A50T ^ 1, where r 
is the scattering time. Under this condition, the peak 
which wc calculated is distinguishable, since the back- 
ground is relatively fiat. However, in order for interband 
peak to dominate over the intraband contribution, a more 
stringent condition must be met. Namely, the products 
A^T, AsoT must exceed {EprY/'^. 

It is worth noting that ideas outlined here have been re- 
cently successfully tested experimentally^ in quasi-one- 
dimensional spin- 1/2 antiferromagnet CS2CUCI4 where 
uniform DM interaction on chain bonds was found to 
split the ESR signal into two peaks already in the 
finite-temperature spin-liquid phase, just as we de- 
scribe above. Many higher-dimensional candidate spin- 
liquid materials, such as two-dimensional k-(BEDT- 
TTF)2Cu2(CN)336 and three-dimensional Na4lr308^, 
have been argued to possess spin-liquid ground states 
with spinon Fermi surface o^'^'^'^i'^^'^^ . DM and other 
kinds of spin-orbit interactions are most certain to be 
present in these materials - this is certainly the case for 
Na4lr308 and highly likely for the organic material as 
well (its close relative, k-(BEDT-TTF)2Cu[N(CN)2]C1, 
does have substantial DM interactioniS) . Our work 
shows that spin-orbit interaction can turn ESR absorp- 
tion experiment into a insightful probe of the unusual 
spinon Fermi surfaces in these and related materials. 
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